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We present transport coefficients (shear viscosity, rj, and bulk viscosity, () for the gluon 
system obtained by the lattice QCD. This is an indispensable calculation towards the 
understanding of "New State of Matter" observed in RHIC. We study the temperature 
regions of RHIC (1.4 < T/Tc < 1.8) and much higher ones up to T/Tc ~ 20. In RHIC 
regions, the ratio of shear viscosity to entropy density, rj/s, is around ~ 0.1 — 0.4, and 
satisfies the KSS bound. At high temperature, rj becomes two or three oder of magnitude 
larger. 

Our calculation has two limitations: (i) the use of the quench approximation, i.e., 
without quark pair creation-annihilation effects on vacuum, and (ii) the use of an ansatz 
for the spectral function. The ffist point has been well studied in calculations of the 
spectroscopy and the phase-transition behavior. To investigate the second point, we 
compare our results with perturbative calculations in high T-regions, and also check the 
effects of the modification of the spectral function on the viscosity. 

1. INTRODUCTION - Matter in deconfinement region 

More than twenty years ago. Gross, Pisarski and Yaffe[|lJ wrote as follows: "Now that 
we possess a theory of the strong interactions, it is natural to explore the properties 
of hadronic matter in unusual environments, in particular at high temperature or high 
baryon density. There are three places where one might look for the effects of high 
temperature and/or large baryon density, (1) the interior of neutron stars, (2) during 
the collision of heavy ions at very high energy per nucleon, and (3) about 10~5 sec after 
the big bang". We are now in an excited era. At RHIC, the confinement /deconfinement 
transition temperature is probably exceeded, and for the ffist time in science history, a 
deconfinement system is created in a laboratory on earth. 

The outcome is very surprising: The matter produced does not look as a quasi-free 
gas as naively expected, but rather is well described as a fluid. In SPS energy regions, 
the hydro- model describes well one-particle distributions, HBT etc., but fails to describe 
the elliptic flow data. This may be not so surprising. Fifty years ago. Landau criticized 
Fermi's statistical model[|2], and noticed 'owing to high density of the particles and to 
strong interaction between them, one cannot really speak of their number' and proposed 
his relativistic hydro-dynamical model[|3j. The ffist quantum field theoretical analysis of 
the applicability conditions of the Landau hydro-dynamical model was reported in Ref . [ 
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In three-dimensional hydro-dynamical calculations to analyze RHIC data, it is assumed 
that the matter produced is a perfect fluid, i.e., its viscosity is zero. This assumption is 
supported by several phenomenological analyses. This also suggests that the new state of 
matter produced at RHIC should be treated as a strongly coupled system: The pertur- 
bative calculation results in 



gHn{fi*/gTy ' ' 

The viscosity is small when g is large. This is understandable because there should be 
sufficient frequent momentum exchange to realize a perfect fluid. Policastro, Son and 
Starinets have shown an example of a strongly coupled theory in which the viscosity is 
indeed very small, i.e., Tj/ s = l/47r[El Ej- They stressed that this is much smaller than 
that of ordinary matter, such as water or liquid helium, and conjectured that this is the 
lowest bount (KSS bound) [[7j. 

It is thus important now to calculate the transport coefficients from QCD, non-perturbatively. 



2. Transport Coefficients on lattice 

On the lattice, the calculation of the transport coefficients is formulated in the frame- 
work of the linear response theory [ |Hl El • 

^ = _ /■ rf3^/ f dtie^(*^-*) r dt'{T^2{x,t)T^2{x',t'))ret. (2) 



Here, {T^uTpa)ret is the retarded Green's function of energy momentum tensors T^^ at a 
given temperature. In the quenched approximation, the energy momentum tensors are 
constructed from only gluonic field strength terms. Bulk viscosity is defined in a similar 
manner. 

Shear viscosity in Eq. Q is also expressed using the spectral function p of the retarded 
Green's function p(a;)[P as 

= vrlim^-^ = vr lim (3) 

It is determined by the shape of the spectral function near a; = 0. 

For evaluating p{uj), we use a well known fact that the spectral function of the retarded 
Green's function at temperature T is the same as that of Matsubara-Green's function. 
Therefore, our target is to calculate Matsubara-Green's function(G^(t„)) on a lattice and 
determine p from it [ E] • 

To determine the spectral function p{uj) from Gp{tn), we adopt the simplest non-trivial 
ansatz, i.e., a Bright- Wigner type ansatz proposed by Karsch and Wyld[|12j, 

PwH = -(7 — ^-7 ■ — — ^) (4) 

As this formula has already 3 parameters, to determine them, the lattice size in tempera- 
ture direction(A^T) should be A^^ > 8. Thus, the minimum lattice size should be 24^ x 8, 
to obtain non trivial results. 
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Figure 1. Shear viscosity as a function of temperature. Its ratio to entropy density in 
RHIC temperature regions (left) and data in physical unit in wide temperature regions 
together with perturbative calculations (right). 



Simulations are carried out using the Iwasaki's improved action and standard Wilson 
action. The simulations are performed at (3 =3.05, 3.3, 4.5 and 5.5 for the improved action 
and at (3=7.5 and 8.5 for Wilson action. With roughly 10^ MC measurements at each (3, 
we determine Matsubara-Green's functions Gpltn)- The errors of G/3 are still large in the 
large t region, however, we fit them with the spectral function p given by Eq.Q. 

The bulk viscosity is equal to zero within the range of error bars, whereas the shear 
viscosity remains finite. 



3. Conclusions 

We may compare our results with the perturbation results of rj in rather high tempera- 
ture regions. In perturbation, bulk viscosity becomes zero[iniEl; whereas shear viscosity 
in the next-to-leading log is given by Eq.(P3). As seen in the right-hand figure of Fig.l, 
in low-T regions, the perturbative calculation becomes inapplicable. At very high tem- 
perature, lattice and perturbative results are satisfactorily consistent with each other. 
Although our result depends on the assumption regarding psw given in Eq. Hi it may be 
a reasonable approximation of dp/ duo at a; = 0. 

Aarts and Resco has proposed an another form of p as [ ITl] 

p{uj)=p'"-{uj)+p'^'^'{uj), (5) 



p'^-^c) = eiuj - h(^/2) + 0.5], (6) 

where dA = — 1 and n{uj) = l/(exp(co'/T) — 1). p^°'^{ijj) is a rational function with 
coefficients as parameter. 

In order to study the effect of p^^^^ on the shear viscosity, 77, we assume that p is given by 
P = Pbw+ P^^^^ ^ where pBw is given by Eq.Q. By changing rrith, the change in r] is studied 
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at P = 3.3 of improved action. When p'^^^'^ is absent ( rrith = oo), 770^=0. 00225(201). If 
mth is set to be 5.0, 3.0 and 2.0, r/a^ becomes 0.00223(0.00191), 0.00194(0.00194) and 

0. 00126.0.00204), respectively. At rrith = 1.8, the contribution from p'^^^h becomes larger 
than Gpitn) of simulation at t = 1, that fit could not be done. Generally, as rrith decreases, 
the contribution from p^^^^ increases and p in the small uj region is suppressed. In this 
case, it results in a decrease in 77. 

We have calculated Matsubara-Green's function and determine the shear viscosity of 
gluon plasma. In the high-temperature region, the agreement of the lattice and perturba- 
tive calculation is satisfactory. The lattice result of r]/ s in T/Tc < 3 is smaller than that 
obtained by the extrapolation of the perturbative calculation and satisfies the KSS bound. 
From the well known relation between the mean free path and viscosity, our results also 
suggest that gluon plasma is strongly interactive. 

Although our results depend on the form of the spectral function psw given by Eq.(jD), 
we think that the quahtative features will not change, because as discussed, our results 
are stable if the high frequency part of the spectral function is included. We think that rj 
and r]/s will not reach 10 times of the present value when more accurate determination 
of the transport coefficients is carried out. 

However, it is important to carry out a more reliable and accurate calculation of trans- 
port coefficients, independent of the assumption regarding the spectral function. To this 
purpose, we are starting the simulation on an anisotropic lattice, to apply maximum 
entropy method. 
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